Devil's staircase in kinetically limited growth 
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The devil's staircase is a term used to describe surface or an equilibrium phase diagram in which various ordered 
facets or phases are infinitely closely packed as a function of some model parameter. A classic example is a 1-D Ising 
model jlj wherein long-range and short range forces compete, and the periodicity of the gaps between minority species 
covers all rational values. In many physical cases, crystal growth proceeds by adding surface layers which have the 
lowest energy, but are then frozen in place. The emerging layered structure is not the thermodynamic ground state, 
but is uniquely defined by the growth kinetics. It is shown that for such a system, the grown structure tends to the 
equilibrium ground state via a devil's staircase traversing an infinity of intermediate phases. It would be extremely 
difficult to deduce the simple growth law based on measurement made on such an grown structure. 



68.55.-a, 81.30.-t ,75.10.Hk 

The original devil's staircase is a footpath between 
Kingshouse to Kinlochleven in Scotland, so called be- 
cause of the huge number of discrete steps between Glen- 
coe and the ridge. In technical usage, the term has been 
used to describe situations in which the number of dis- 
crete steps within a finite range becomes formally infinite. 
Examples include inter alia the formation of facets of a 
crystal antiferroelectric, smectic and lyotropic liq- 

uid crystals , magnetic structure in cerium monop- 

nictides |(| and granular media Usually the stair- 
case emerges from the interplay between long-range re- 
pulsive (antiferromagnetic) and short range attractive 
(ferromagnetic) forces, with transitions between stable 
phases appearing as the relative strengths of the interac- 
tions are altered. The precise form of the interactions is 
not crucial [||. 

The drive toward nanofabrication has led to a tremen- 
dous interest in growing multilayer structures. In typical 
methods such as molecular beam epitaxy or chemical va- 
por deposition careful control of the composition of the 
deposited material is required to create complex artificial 
structures. Without such careful control non-periodic 
structures tend to form. By contrast, some structures 
of technological interest such as quantum dots may self- 
assemble, and understanding the local equilibria which 
govern growth is crucial. In this letter layer-by-layer sur- 
face growth for a simple model is shown to yield a devil's 
staircase structure. This suggests that for a wide class 
of systems the expected structure grown at zero temper- 
ature is aperiodic, and might easily be misinterpreted as 
disordered. This apparent disorder is not real, but arises 
from the locally stable structure being dependent on the 
thickness of the film, and there being an infinite number 
of locally stable structures. 

Specifically, the Hamiltonian for our model is equiva- 
lent that considered by Bak and Bruinsma ifflj 
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This model describes a situation in which each layer 
can be one of two types cr, = ±1. The first term gives 
Cj = — 1 a lower formation energy than <Xj = +1, while 
the second term gives a long ranged repulsion between 
like-layers. 

Previous work has concentrated on the devil's stair- 
case as an equilibrium phenomenon, and searched for the 
thermodynamic ground state. Here, by contrast, the dy- 
namics of growth are considered, spins being added to 
the system so as to minimise the energy, but then being 
fixed forever as further layers grow. 

Many physical systems can be mapped onto this 
Hamiltonian: a simple example is a line of charges in 
an external field. The same Hamiltonian describes a sit- 
uation where the Oj represent the separations between 
layers rather than the layers themselves. Now the first 
term indicates that it costs less energy to grow either type 
on a similar layer, while the second term again indicates 
long-range repulsion(attraction) between likc(unlikc) lay- 
ers. This might describe a system where epitaxial growth 
was favored, but generates a long range strain field which 
needs to be periodically relieved. 

Alternately it may describe a situation such as silicon 
carbide growth or stacking of close-packed planes, 

where each layer is locally either ABA or ABC stacked 
depending on its neighbors. Now a represents the relative 
orientation of adjacent layers. In close packed layer (AB) 
and interlayer (er) notation equivalent stacking sequences 
for fee containing a growth fault are 



ABCABCACBACBA. 

. + + + + + + - + + + + + + . 



(2) 
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Notice that a single fault of this type cannot be accom- 
modated within periodic boundary conditions. This led 



1 



Bak and Bruinsma to postulate that the actual defects in 
the devil's staircase are fractional, since more than one 
must be created together. In the growth case there is 
no such constraint: this is equivalent to the difference 
between intrinsic and extrinsic stacking faults in close 
packed materials (which can arise from removal or inser- 
tion of a plane) and growth defects (basal plane twins) 
which reverse the sense of stacking and can be generated 
only by finite shear of the entire sample or during growth. 
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Finally, the model can describe a simple history- 
dependent system, where the state of the system depends 
on a sum over its historical values. In this case the 
"layer number" should be interpreted as a time rather 
than space dimension. 
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For the growth dynamics, one simply considers adding 
the n+lth. layer to the preexisting n layers with whichever 
spin reduces the energy. This can be determined entirely 
by the sign of the local potential. 
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In zero temperature case considered here, if V n -\-i is 
positive, the next added layer a n +i = — 1 5 otherwise 
(j rl+ i = +1. furthermore, the final structure is uniquely 
defined and while it may appear random, it has zero en- 
tropy. 

At thermodynamic equilibrium, or asymptotically for 
the growth dynamics, this model (equation 1) has two 
simple limiting cases. For weak long-range interactions 
defined by 



FIG. 1. Self-similar devil's staircase of phases reached 
asymptotically in growth with v — 2. Plotted are the value of 
A and the mean values of a evaluated over the final 2520 lay- 
ers of a 300000 layer sample. 2520 (= 2x3x2x5x7x2x3) is 
chosen because it is commensurate with all periodicities from 
one to ten layers, and with 12,14,15,18,20,21 etc). For these 
case the plotted value of < a > is exact, for others it may be 
out by up to 0.04%. 



In the case of growth, we find that a convenient pa- 
rameter to monitor is < V > m , the rolling mean value 
of the potential between V n - m arL d V n wherever a layer 
type a n = +1 is grown. For the case of v = 2 the asymp- 
totic value of < V >2520 plotted against A picks out the 
conventional devil's staircase behavior (Figure 1). 
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<7„ = —1 for all n, meanwhile for small A alternating 
behavior a n = (—1)™ is observed. For intermediate val- 
ues of A, the devil's staircase of phases is recovered in 
the asymptotic limit (figure 1) jjj. 



Our interest lies in the convergence of the structure 
with layer number - physically how thick a film must 
grow to recover bulk behavior. Again this can be mon- 
itored using < V > m , now plotted against layer num- 
ber. For v = 2 the growth converges fairly rapidly to 
the equilibrium value, the effective screening of the sur- 
face is fast compared with the integer layer spacing. For 
smaller v convergence is slow: Figure ^| shows the case 
of v = 1,A — 1: now the screening is sufficiently slow 
that a wide range of different phases from the "devil's 
staircase" are actually observed over a number of layers. 
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FIG. 2. Mean surface potential averaged over the preced- 
ing 72 layers before growing the i th a — +1 layer. Highly 
ordered regions correspond to short-period phases which are 
stable over a significant range of thickness. The slope of the 
graph shows that even within these regions, equilibrium has 
not been reached and ultimately the order breaks down as a 
new phase is stabilized. A long range trend towards an asymp- 
totic value of < V >72 is observed. This figure generated for 
v = l, A=l, V(0)=0 



3 



-0.56 



-0.57 



A 
> 
V 



i ! I 
i I 



I i i 
i i 



■ \ 

\\ \ 

A \\ \ - 

I i\\ Vi 
1 ' SB \\\ « ■ 



\\ 



w; 



-0.58 



-0.59 



20000 



40000 



60000 80000 
Number of layers 



100000 



120000 



FIG. 3. Detail from figure 2, illustrating the ordered nature of of the potential variation for a region where long-period 
phases are stable. The long-period repeats give rise to multiple values of < V >72 within the same phase, hence the multiple 
branches. This multiplicity can be reduced or eliminated by considering < V >2520 at the expense of smearing out details. In 
this respect the figure is not self-similar. The (arbitrarily-chosen) 72-layer averaging is significant on the scale of this figure, 
manifesting itself for the phases which are stable over more than 72 layers as an initial increase and subsequent curvature in 
< V >72- The pure phase behavior is typically a linear decrease of < V >j2 with N, as seen after 72 layers for those phases 
which are stable over a sufficiently long period. 



The actual phases and "screening" effect are illustrated 
by the running average of the ratio of<r = +lto<7 = — 1 
over the previous 72 layers (72 is chosen such that phases 
with period 2,3,4,6,8,9 etc will give a constant value). To 
further reduce oscillation, the ratio is printed out only at 
layers with a = +1. Substantial single phase regions can 
be seen, together with shorter transitional regions. The 
overall trend towards a limiting value can be seen. 

Each phase is stable only over a finite number of layers. 
Since the range of stability is inversely proportional to the 
period of the structure the thickness over which some 
long-period structures are stable will be shorter than 
the periodic repeat distance of these structures. Con- 
sequently, they cannot be unambiguously identified. 

The long term trend of Figure 2 is toward < V > m = 
—0.5. This corresponds to equal numbers of a = ±1 



which give a mean field value of V = averaged over 
all layers, and < V > m = —A/ 2 = —0.5 averaged over 
the (7 = +1 layers only. A very curious phenomenon ob- 
served in Figure 2 is that the ordered phases show antis- 
creening behavior: the mean value of < V > m moves 
away from the asymptotic value for most of range of 
the ordered phase. Thus while the overall trend is an 
asymptotic increase of < V > m , within any given phase 
d < V > m /dN is negative [l^ . Evolution towards the 
asymptote is achieved via the boundaries between the 
phases, rather than the screening by the phases them- 
selves. 
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FIG. 4. Mean value of <x for v = 1, A = 1. The upper 
line is averaged over the preceding 2520 layers, picking out as 
straight lines phases of repeat period as in figure 1. Transi- 
tions between these short period phases are characterized by 
longer period phases, which are shown by thick black lines (ac- 
tually representing a rapid fluctuation between K/2520 and 
(K+l)/2520 for integer K). The lower line is averaged over 
the whole layer, and shows the very slow monotonic growth 
of < a > towards the asymptotic value of < a >= 0). The 
increase in the mean value of 1/ < a > is logarithmic, a 
reasonable fit to the graph being < a >= [7 — 2. 9 In TV] _1 . 
For higher values of A the convergence of < a >^ is even 
slower. For v > 1 the asymptotic value of < a > is non-zero: 
it is a phase from the devil's staircase dependent upon A. 

In mapping onto a real growth process, the value of 
A is determined by the material being deposited, how- 
ever it may be possible by choice of substrate or external 
applied field to control the initial value of the potential. 
By doing so, the density profile of a = +1 may be var- 
ied. This is not straightforward however: if the initial 
condition is compatible with the equilibrium structure a 
perfect multilayer can be grown, if not the concentra- 
tion will traverse all possible phases with a = +1 density 
intermediate between the starting and equilibrium ones. 
There are an infinity of these phases, but their thick- 
ness must take an integer value, thus not all phases can 
actually be observed. If the interlayer spacing is taken 
to be of atomic dimension, a perfectly grown (i.e. zero- 
T, zero entropy) film of even a millimeter thickness may 
not reach equilibrium and will appear disordered to any 
experimental probe. 

By interpreting the layer number N as a time rather 
than a thickness, this type of growth-kinetic model also 
provides a simple model of history-dependence. In many 
social phenomena, decisions are made for one of two 
courses of action based on the evidence of past behavior 



with more recent evidence having a stronger weight [[F3| . 
Here the model already contains enough complexity to 
behave counterintuitively: the long term trend of increas- 
ing < V > m is opposite in sign to the d < V > m /dN 
measured over the stable phases, despite the fact that 
formally the devil's staircase provides a stable phase at 
all N, and by implication d < V > m /dN is negative ev- 
erywhere. Of course, the spin-Ising model is a gross over- 
simplification of any real decision making process: this 
only serves to emphasize the non-trivial relation between 
the model and its behavior, and the difficulty for mea- 
surement when d < V > m /dN is negative everywhere 
while AV/ AN is positive. 

In sum, the growth kinetics of a simple model exhibit- 
ing long range antiferromagnetic and short range fer- 
romagnetic ordering have been studied. This has been 
shown to exhibit logarithmically slow equilibration to the 
equilibrium structure, passing through a formally infinite 
number of intermediate phases which form a devil's stair- 
case. 
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